FERMILAB-Conf-04/096-E 



B Physics at D0 



Jan Stark 
on behalf of the D0 Collaboration 



Laboratoire de Physique Subatomique et de Cosmologie, Grenoble, France 



The Fermilab Tevatron (pp), operating at y/s = 1.96 TeV, is a rich source of B hadrons. 
The large acceptance in terms of rapidity and transverse momentum of the charged particle 
tracking system and the muon system make the upgraded Run II D0 detector an excellent 
tool for B physics. In this article, we report on selected physics results based on the first 
250 pb _1 of Run II data. This includes results on the X(3872) state, semileptonic B decays, 
B hadron lifetimes, flavour oscillations, and the rare decay B s — > 



1 Introduction 



The D0 Collaboration has implemented a major upgrade of their detector for Run II of the Fer- 
milab Tevatron. Certain components of this upgrade have significantly improved DCs B physics 
potential compared to Run I. The upgrade is described in more detail in Ref. [1]. One of the 
aspects most relevant for B physics is the new charged particle tracking system located in a 
2 T field produced by a superconducting solenoidal magnet. The tracking system consists of a 
scintillating fibre tracker and an inner tracker based on silicon detectors. It provides coverage 
over a wide range in pseudo-rapidity (\rj\ < 2, with potential for extension up to \rj\ < 3) and 
allows the efficient reconstruction of tracks, including down to low transverse momentum of 
Pt = 180 MeV/c. Due to the limited radius of the fibre tracker (imposed by the geometry of 
calorimeter, reused from Run I), the track pT-resolution is not as good as, e.g., in CDF, which 
translates into less good resolutions on the mass of fully reconstructed B hadrons. Typical 
mass resolutions for decay modes^ like B + — > J/ip(— * p, + fi^)K + are 40 MeV/c 2 , which still 
allows a clean identification of B hadrons. Track impact parameters do, in the plane trans- 
verse to the beam axis, with respect to the primary vertex are measured with high precision 
(c(drj) = 50 fim for tracks with j>t — 1 GeV/c, improving asymptotically to a (do) = 15 ^m for 



"Throughout this article, references to a hadron or to a decay reaction also imply their charge conjugate. 



tracks with pr > 10 GeV/c), which allows precise measurements of the time-dependence of B de- 
cays (lifetimes, flavour oscillations, CP violation, etc). Another key ingredient of the upgraded 
D0 detector is the muon detector and its associated trigger system. The muon detector, located 
outside the calorimeter, relies on three layers of scintillation counters and drift chambers, one 
of which is placed in front of toroidal magnets, the other two behind the toroids. Muons with 
momentum greater than ~ 1.4 GeV/c (~ 3.5 GeV/c) penetrate the calorimeter (toroid magnet) 
at rj = 0. The muon system provides wide coverage over \rj\ < 2.0, including at the trigger level. 
The data used for the measurements reported here have been recorded with single and di-muon 
triggers. The single muon triggers are prescaled away at the highest instantaneous luminosities. 
Most of the results reported in this article are based on 250 pb _1 of data (before prescales) 
collected during the period from April 2002 to January 2004. Some of the results are based on 
earlier subsets of this data sample. 

2 Production and decay characteristics of the X(3872) 

The Belle Collaboration has recently discovered a new particle, the X(3872), with a mass of 
3872.0 ± 0.6 (stat) ± 0.5 (syst) MeV/c 2 , produced in B ± -» X(3872)K ± , X(3872) -» J/^vr+vr- 
decays [2]. The existence of the X{2>872) state has been confirmed, also using the J/ipTr + ir~ 
mode, in pp collisions by CDF [3]. At this time, it is still unclear whether this particle is a 
cc state or a more complex object; see, e.g., Refs. [4-6]. 

We perform a more detailed study [7] of the production and decay characteristics of the X(3872) 
in pp collisions. The cc state tp(2S) has a mass m(ij>(2S)) = 3685.96 ± 0.09 MeV/c 2 [8] close to 
the X(3872), and has the same decay mode J/ilm + ir~ . It therefore provides a good benchmark 
for comparison with the X(3872). Using 230 pb" 1 of data, we reconstruct 522 ± 100 A(3872) -> 
J/if)TT + ir~ decays. We split our samples of ip(2S) and A(3872) decays into two subsamples 
according to transverse momentum p^(J/tp7r + -K~), and compare the event-yield fractions f^s) 
and .fx (3872) w hh 

_ N(P) at pr(j/ijm+ir-) > 15 GeV/c 2 
~ N(P) at pt(J/^tt + tt-) < 15 GeV/c 2 ' 

where P = 4>{2S) or X(3872). The result is shown in Fig. 1. We repeat this study for splittings 
according to the rapidity y of the J/tpir + TT~ system, the helicity of the tt + tt~ (/i + /i~) system, 
the effective proper decay length and isolation (Fig. 1). The helicity of the 7r + 7r~ system 
is inferred by boosting one of the pions (muons) and the J/tpir + ir~ system into the di-pion (di- 
muon) rest frame and measuring the angle 9 n (6^) between them. The effective proper decay 
length dl is defined as the distance in the transverse plane from the primary vertex to the decay 
vertex of the J/Y>, scaled by the mass of the J/^ir + Tr~ system divided by Pt(J /^tt + it~). The 
isolation is defined as the ratio of the J/?/>7r + 7r~ momentum to the sum of the momentum of 
the J/Vw + 7r~ and the momenta of all other reconstructed charged particles within a cone of 
radius AR = 0.5 in rj/(f> space about the direction of the J/^tt + tt^ momentum. In summary, 
no significant differences between the X(3872) and the cc state i){2S) are found when the data 
are separated according to any of these production and decay variables. 

3 Precision measurement of t(B + )/t(B°) 

The spectator quark model predicts that the two charge states with one heavy quark Q (Qu and 
Qd) have the same lifetime. Deviations from this simple picture have been evaluated from first 
principles of QCD using Operator Product Expansion (an overview is presented, e.g., in Ref. [9] 
and references therein) and are expected to be proportional to l/[m(Q)] 3 . Using inputs from 
quenched lattice QCD, recent NLO calculations [10] yield t{B+)/t(B°) = 1.053 ±0.016 ±0.017, 
but unquenching effects could be sizeable. 
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Figure 1: Comparison of event-yield fractions for X(3872) and ^){2S) in the regions: (a) Pt(J '/ipir^ tt ) > 
15 GeV/c; (b) \y(J/ipir + ir~)\ < 1; (c) cos^) < 0.4; (d) effective proper decay length, dl < 0.01 cm; (e) isola- 
tion = 1; (f) cos(0 M ) < 0.4 . 



We measure t(B + )/t(B°) in 250 pb^ 1 of data using large samples of semileptonic B + and 
B° decays collected using single muon triggers. Semileptonic B decays are reconstructed in the 
final state D°(—> K + ir~)fi + X, where X denotes any additional decay products. Figure 2 shows 
the m{K + Ti~) spectrum for K + ir~ /i + events after pt and tracking hit requirements on the tracks 
and cuts on the significance of the impact parameter of the K + and tt" tracks, the significances 
of the distances in the transverse plane between the primary vertex and the D° and B vertices, 
the angle in the transverse plane between the D°[i + momentum and the direction of the D°fx + 
vertex w.r.t. the primary vertex, as well as the requirement 2.3 < m(L>° / u + ) < 5.2 GeV/c 2 . 
The signal in the D° peak contains 109/c events. The next step is to identify the subsample 
of events in which the D° originates from a D*~ — > D°tt~ decay. Due to the small D*~ / D° 
mass difference, the transverse momentum of the pion tends to be low, and for small transverse 
momentum of the D*~, the slow pion reconstruction efficiency is limited by tracking acceptance. 
For each candidate, an additional pion with charge opposite to the charge of the muon is 

searched for. The mass difference Am = m(D ir~) — m(D°) for events with such an additional 
slow pion is shown in Fig. 3. The peak corresponds to D*~fi + events. 

The slow pion reconstruction allows the sample to be split into a "D° sample" (D°(— > K + it~)[i + X 
events in which no slow pion is found) and a "D* sample" (D°(— ► K + tt~)tt~ /j, + X). The D° sam- 
ple is enriched in B + decays, while the D* sample is enriched in B°. This separation be- 
tween the B species gives us access to the lifetime ratio t(B + )/t(B°). We further require 
Pt(D°) > 5 GeV/c. This cut reduces the event yields by ~ 35 %, but it also ensures a good slow 
pion reconstruction efficiency (88 %, independent, within uncertainty, of the B meson decay 
length) and therefore a good separation between the D° and D* samples. The relative abun- 
dances of different species of B hadrons in the D° and D* samples are evaluated from measured 
branching fractions and their uncertainties, and the reconstruction efficiencies for the different 
decay chains. Isospin relations are used to infer some unmeasured branching fractions from 
measured ones. The D° sample contains 82 % B + , 16 % B° and 2 % B s , while the D* sample 
contains 86 % 5°, 12 % B + and 2 % B s . 

As the semileptonic B decays in our samples are only partially reconstructed, the proper decay 
time of a B candidate can only be inferred approximately from the measured decay length in 
the transverse plane It- For each event, we calculate the visible proper decay length (VPDL) 
defined as x m = It ■ - , and we require a(x m ) < 200 /xm. The expected VPDL distribu- 
tion for a given B species and a given decay chain is obtained from the expected proper decay 
time distribution and the distribution of K = pt(D° fx) /pr(B) ( "fT-factor" , determined from 
simulated events). The lifetime ratio t(B + ) /t{B ) could be extracted from fits to the VPDL 
distributions of the D° and D* samples. We chose a different technique. We group the events 



into 8 bins of VPDL, and measure the ratio 



ri a^V+) 

separately in each bin i. The event yields Ni(D*~ /i + ) and Ni(D°fi + ) are determined from 
separate fits to the m(K + ir~) mass spectra (similar to the one in Fig. 2) of the D° and D* 
samples in each bin. The D* sample contains a small contribution of combinatorial background 
from events with a true D°, which is counted as signal in the m(if + 7r~) fit. This contribution is 
subtracted using the m(K + 7r~) distribution of wrong charge combinations D°tt + in each VPDL 
bin and a correction factor C to take into account differences in the combinatorial background 
between right and wrong charge combinations. The correction factor C = 1.22 ± 0.04 is derived 
from the full sample. The lifetime ratio t(B + )/t(B°) is determined from a fit to the 8 values rj. 
An advantage of this technique compared to fits to VPDL distributions is that the background 
is properly subtracted in all VPDL regions and no parameterisation of the VPDL distribution 
of background is needed in the fit. This technique also provides an elegant way to take into 
account the decay length significance cuts described above. Additional inputs to the fit are the 
relative B contents of the D° and D* samples, taking into account the effect of reconstruction 
efficiencies (in particular the slow pion reconstruction efficiency) , the distributions of K-factors 
for different decay chains, the decay length resolution (from simulation) and the lifetime t(B + ) 
which is fixed to the PDG value t(B+) = 1.674 ± 0.018 ps [8]. 

Figure 4 shows the measured values of the as well as the result of the fit. We obtain the 
preliminary measurement 

t(B + )/t(B°) = 1.093 ± 0.021 (stat) ± 0.022 (syst) . 

The systematic uncertainties are dominated by the uncertainty on a possible decay-time de- 
pendence of the slow pion reconstruction efficiency, uncertainties in the relative reconstruction 
efficiencies for different B hadron decay chains, the uncertainty in Br(i? + — > D*~ir + fi + isX) 
which is used in the estimate of the sample compositions, the if-factors from simulation and 
the uncertainty in the decay length resolution differences between the D° and D* samples. This 
preliminary D0 measurement is one of the most precise measurements of t(B + )/t(B°), and it 
is consistent with previous measurements [8]. It is also consistent with the most recent prelim- 
inary CDF measurement, t(B + )/t(B°) = 1.080 ± 0.042, based on fully reconstructed hadronic 
B -» J/tj)X decays in 240 pb" 1 [11]. 

4 B° I B° mixing with opposite-side muon tagging 

We use our sample of semileptonic B° — > D*~ fi + u decays to measure the B°/B° mixing fre- 
quency Amrf. Compared to the lifetime measurement, we need an additional key ingredient: the 
ability to identify ("tag") the initial flavour of the B meson. Opposite-side muon tagging relies 
on the charge of muons among the decay products from the "other b" from bb pair production to 
infer the initial flavour of the reconstructed B. Typically, the bb pair leads to two-jet events with 
two B mesons in two back-to-back jets. The final state can also contain a b baryon and/or more 
than two jets. In addition, if the two 6-jets originate from the flavour excitation or gluon splitting 
processes, the angle between the two jets is not necessarily close to 180 degrees, but varies within 
a rather wide range. We select tagging muons with pr > 2.5 GeV/c and cos(A0) < 0.5, where 
Acj) is the difference between the (f) directions of the momentum of the D°n + system from the 
reconstructed B and the momentum of the tagging muon. This algorithm provides an estimate 
of the initial flavour for (4.76 ± 0.19) % of all reconstructed B°/B° mesons. 
The mixing frequency Am^ has been measured with high precision by many experiments, most 
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Figure 2: Invariant mass m(K + Tv~) of the 
7r system in K it fi events. The solid lines 
represent the result of a fit to the sum of a gaus- 
sian distribution to model the D° — > K + ir~ peak 
and a polynomial background model. Data points 
in the low mass peak, corresponding to partially re- 
constructed D° — > K + tt~ X decays, are assigned a 
zero weight in this fit. 
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Figure 3: Mass difference m(D°7v^) — m(D°) for 
D°n-fi+ events with 1.80 < m(D°) < 1.90 GeV/c 2 . 



recently at the U-factories. We use the study of B° / B° mixing to benchmark our flavour tagging. 
A fit to the observed asymmetry between the number of mixed and unmixed events as a function 
of VPDL allows us a simultaneous measurement of Am^ and the tagging purity k = jj^-+jf^, 
where + iVw is the number of events for which the tagging algorithm has provided an es- 
timate of the initial flavor, and Nn is the number of events in which the estimate was correct. 
This allows us to measure the performance of our flavour tagging in data before we use it in 
more difficult studies of, e.g., B s /B s mixing and CP violation. The measurement of Am^ also 
constrains more exotic models of b production in hadron collisions (see, e.g., Ref. [12]). 
Using 250 pb -1 of data, we obtain the VPDL-dependent asymmetry between mixed and unmixed 
events shown in Fig. 5 and the preliminary measurements 

Am d = 0.506 ± 0.055 (stat) ± 0.049 (syst) ps" 1 and k = 73.0 ± 2.1 (stat) ± 0.8 (syst) % . 

The systematic uncertainties on Amd are dominated by the uncertainties in the method used 
to fit the number of signal events in each VPDL bin. The systematic uncertainties on k are 
dominated by the uncertainty in the B s contribution to the signal sample, on the proper decay 
time resolution and in the alignment of the tracking system. The measured value of Amd is 
consistent with the world average [8]. 

This kind of analysis is currently being repeated for other, complementary, tagging algorithms 
based on opposite-side electrons, jet charge and same-side tags. These additional taggers will 
further improve our overall flavour tagging efficiency. 



5 Observation of semileptonic B decays to narrow D** mesons 

Orbitally excited states of the D meson are usually referred to as D** mesons. In the simplest 
case, D** mesons consist of a charm quark and a light quark in a state with orbital angular 
momentum L = 1. In the limit m(c) 3> Aqcd ; the spin of the charm quark decouples from the 
other angular momenta, and the angular momentum sum j = S + L of the light quark spin S 
and the orbital angular momentum L is conserved. For L = 1, one therefore expects one doublet 
of states with j = 3/2 and another doublet with j = 1/2. The two j = 3/2 states are denoted 
D® and D^ . In contrast to the j = 1/2 states, the and decay through a D-wave and 
are thus expected to be narrow (20 to 30 MeV/c 2 , while the j = 1/2 states decay through an 
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Figure 4: Measured values of n 
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the 8 bins of visible proper decay length (VPDL). 
The solid histogram represents the result of the fit 
to the data. 




Figure 5: Asymmetry between the number of mixed 
(iV osc ) and the number of unmixed (JV non -° sc ) 
events as a function of visible proper decay length 
(VPDL). The solid line represents the result of the 
fit that is used to measure Am<j and the tagging 
purity k. 



S-wave and are expected to be several hundred MeV/c 2 wide). For a recent theoretical overview 
of the topic, see, e.g., Ref. [13] and references therein. 

The Di and D^ have been observed and studied in several experiments, most recently by 
Belle [14] and BABAR [15] in the hadronic decay B — > D**°ir~. We study D^, produced in 
semileptonic B decays. Our measurements further constrain D meson spectroscopy and improve 
the understanding of the samples of semileptonic B decays we use as a tool for studies of the 
time-dependence of B decays (lifetimes, mixing, etc). Figure 6 shows the m(D*ir) spectrum 
we observe in candidate B — > (D?, D^fM ~vX — > D* + ir~ \i~vX decays using 250 pb" 1 of data. 
The excess in right-sign combinations can be interpreted as combined effect of the D\ and 
From topological analyses at LEP [8], we know the inclusive rate 

Bt(B -> D* + -K-fi-vX) = 0.48 ± 0.10% . 

Our preliminary result constrains the resonant contribution 

Bi(B {D%,D?}/TvX) • Br({L>?, Df} -» D*+tt-) = 0.280 ± 0.021 (stat) ± 0.088 (syst) % . 

The branching ratio for B decays to D** has been determined using a normalisation to the 
decays B — > D* + (— > D°ir + )fi~h i X, whose branching ratios are well known [8]. The systematic 
uncertainties are dominated by the uncertainty on the pt dependence of the soft pion recon- 
struction efficiency, as well as the effects of possible contributions from wide resonances and 
the effect of possible interference effects between the D® and D^ on the tails of the D** mass 
distribution, which are not taken into account in this preliminary analysis. 



6 Fully reconstructed B decays 

D0's good tracking acceptance and efficiency, including for tracks from, e.g., — > tt + tt^ or A — > 
p + n~ that originate far away from the primary vertex, also allow the exclusive reconstruction of 
hadronic B decays. The example of B° — > J/f/j(— »• fi + p~)Ko i is shown in Fig. 8. Using this decay 
mode in 114 pb -1 of data, we obtain the preliminary measurement t(B°) = 1.56lg25 (stat) ± 
0.13 (syst) ps, consistent with the world average [8]. Figure 8 also shows the A& — ► J/ip(— > 
/i + /i~)A(^ p + Ti~) signal. A lifetime measurement in this channel is currently being finalised. 
Another signal with rich physics potential, B s — > J/ifj(— > n + fj>~)<p(— > K + K~), is also included 
in Fig. 8. 
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Figure 6: Invariant mass m(D*n) for candidate B — > 
(Di,L>|°)/i _ !/X -> D* + tt~ fi~vX events. The points 
correspond to D*7r combinations with opposite charges, 
and the solid histogram corresponds to same charge com- 
binations. The total number of events in the peak is 
523 ± 40, estimated as difference in the number of events 
for right charge and wrong charge combinations with 
2.35 < m(D*7r) < 2.55 GeV/c 2 . 
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Figure 7: Definition of the signal 
and sideband regions in the B s — > 
fi + fi~ sensitivity analysis, and remain- 
ing background in 180 pb _1 of data af- 
ter preliminary cut optimisation. 



7 Sensitivity analysis of the rare decay B s — ► M A* 

As discussed, e.g., in Refs. [11,16,17] and references therein, the rare decay B s — > fi + n~ is a 
promising window on possible physics beyond the Standard Model. We study the sensitivity 
of D0 to this decay using 180 pb" 1 of data. The main discriminating variables used to suppress 
background are the isolation of the muon pair, the transverse decay length \l xy \ of the ver- 
tex w.r.t. the primary vertex, and the opening angle between I and p(fi + The fJ> + [i~ mass 
spectrum after cuts is shown in Fig. 7. Based on the number of events observed in the mass 
sidebands (Fig. 7), we expect 7.3 ± 1.8 background events in the signal region (Fig. 7), and, in 
the absence of an excess over the expected background, a limit of Bv(B s — > < 1.0 • 10~ 6 

at 95 % CL. This expected limit takes into account both statistical and systematic uncertain- 
ties. The expected signal has been normalised to an abundant sample of fully reconstructed 
B — > J/ipK^ decays. Possible contributions from B& — > jU + //~ are expected to be small com- 
pared to B s — > and are neglected here. The signal box has not been opened yet, i.e. the 
number of events in the signal region has not been looked at. The event selection is currently 
being reoptimised, and a very competitive limit on a slightly larger data sample is expected on 
a short timescale. 



8 Conclusion and outlook 

The large acceptance in terms of rapidity and transverse momentum of the charged particle 
tracking system and the muon system, including the muon trigger, make the upgraded Run II 
D0 detector an excellent tool for B physics. In this article we have presented selected physics 
results based on up to 250 pb _1 of Run II data. Large samples of semileptonic B decays have 
been accumulated using single muon triggers. The excellent impact parameter resolution of the 
tracking system allows precise measurements of the decay-time dependence of B decays, and 
one of the world's most precise measurements of t(B + )/t(B°) has been performed. The muon 
reconstruction capabilities also allow pure flavour tagging, which is another key ingredient for 
the study of flavour oscillations and CP violation. The good tracking efficiency down to very low 
transverse momenta has made possible the observation of semileptonic B decays to D** mesons. 




Figure 8: Fully reconstructed b hadrons in three selected decay modes. 



Large yields of hadronic decays of b hadrons to J ftp — > are obtained using di-muon 

triggers. Di-muon trigger data have also been used for a detailed study of the characteristics 
of the A (3872) state and a search for the rare decay B s — > fi + . Other hadronic decays 
are collected by triggering on single opposite-side muons. Events recorded in this way are 
particularly powerful for measurements that require flavour tagging, as the tag is provided by 
the muon that fired the trigger. In addition, the proper decay time distribution of these events is 
not biased by the trigger. Such events will be included in future analyses. Other improvements 
are expected from the future use of dE/dx information from the silicon tracker, which provides 
kaon and proton identification at very low momenta. In combination with information from 
the tracking system, D0's calorimeter is starting to be used to identify non-isolated low-px 
electrons. Together with the calorimeter, the preshowers (which provide finer granularity than 
the calorimeter) are used to identify 7r° decays. In summary, the exciting times of the competitive 
D0 Run II B physics program have started. 
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